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ABSTRACT 

The two brightest elliptical galaxies in the Eridanus A group, NGC 1400 and NGC 
1407, have been observed in both the Washington T\ and Kron-Cousins / filters to 
obtain photometry of their globular cluster systems (GCSs). This group of galaxies is 
of particular interest due to its exceptionally high M/L value, previously estimated at 
~ 3000/;,, making this cluster highly dark-matter-dominated (Gould 1993). NGC 1400's 
radial velocity (549 km/s) is extremely low compared to that of the central galaxy of 
Eridanus A (NGC 1407 with v@ = 1766 km/s) and the other members of the system, 
suggesting that it is a foreground galaxy projected by chance onto the cluster. Using 
the shapes of the globular cluster luminosity functions, however, we derive distances of 
17.6 ± 3.1 Mpc to NGC 1407 and 25.4 ± 7.0 Mpc to NGC 1400. These results support 
earlier conclusions that NGC 1400 is at the distance of Eridanus A and therefore has a 
large peculiar velocity. Specific frequencies are also derived for these galaxies, yielding 
values of S N = 4.0 ± 1.3 for NGC 1407 and S N = 5.2 ± 2.0 for NGC 1400. In this and 
other respects, these two galaxies have GCSs which are consistent with those observed 
in other galaxies. 

Subject headings: Globular cluster systems, galaxy clusters, dark matter 
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1. Introduction 

NGC 1400, a bright E0/S0 galaxy located in a large 
sub-condensation of the Eridanus cluster, has recently 
drawn attention to itself and its curious environment. 
The southern Eridanus group of galaxies extends between 
3 h 15 m £ R.A. £ 4' 1 and -26° £ Dec. £ -15°, and 
exhibits a large degree of internal sub-clustering (see Fig- 
ure 1 from Willmer et al. 1989). Its most concentrated 
clump of galaxies is Eridanus A (3 h 40 m , — 19°), which in- 
cludes roughly 50 galaxies within a radius of ~ 0?8. 

At an estimated distance of 16.4 Mpc (Tonry 1991), 
the large central E0 galaxy of the Eridanus A sub-cluster 
(NGC 1407) has a heliocentric velocity typical of the ma- 
jority of the members of the group. NGC 1400, however, 
reveals an anomalously low redshift (549 km/s) compared 
with those of its neighbouring galaxies, causing some un- 
certainty as to its true distance. Table gj lists the members 
of Eridanus A for which heliocentric velocities have been 
determined. The mean radial velocity of these galaxies is 
1643 ± 95 km/s, or 1765 ± 93 km/s if NGC 1400 is not 
considered a member. 

This straightforward question as to the true distance of 
a potential Eridanus A member has resulted in what may 
prove to be an extraordinary find: the sub-cluster has 
an abnormally high M/L value. Of the 10 Eridanus A 
galaxies listed in Table |l|, the two early-type galaxies 
NGC 1407 and NGC 1400 account for nearly 80% of the 
light. Adding the blue luminosities of the Eridanus gal- 
axies within 0?8 of NGC 1407 (including NGC 1400) and 
calculating the mean of four virial mass estimators, Gould 
(1993) finds M/L ~ 3125ft. 

This value is anomalously large; typically, groups or 
clusters of galaxies have M/L values of several hundred, 
making Eridanus A one of the darkest clusters known. 
Considering only the virial theorem mass, the exclusion 
of NGC 1400 from the cluster mass calculation reduces 
M/L by roughly a factor of 2. Consequently, even if NGC 
1400 has not virialized within the cluster, Eridanus A's 
M/L is still excessively high. 

If NGC 1400 is at the distance of Eridanus A, it may 
yet have originated from some other location, giving rise 
to a peculiar motion which is not due to the mass con- 
centration of the sub-cluster. Since there are no other 
nearby mass concentrations large enough to account for 
NGC 1400's motion, this possibility is unlikely (see Gould 
1993). Therefore, if NGC 1400 is at the distance of Eri- 
danus A, it is likely to be dynamically associated with the 
sub-cluster. Is NGC 1400 indeed at the distance of Eri- 
danus A, or is it instead a foreground galaxy mistakenly 
associated with the group due to projection effects? 

Five independent studies have attempted to determine 
whether NGC 1400 is at the distance of Eridanus A. The 
results of four of these studies are summarized in Table 2 



(see discussion in Gould 1993 and references therein) . One 
additional piece of evidence presented by Gould is the 
globular cluster luminosity function (GCLF), credited as 
showing that NGC 1400 and NGC 1407 are at the same 
distance. However, these GCLFs are secondary results 
obtained by Tonry (1991) during the surface-brightness- 
fluctuation (SBF) analysis, and no rigorous investigation 
of the globular cluster systems of the two bright galaxies 
in Eridanus A has yet been performed. 

We therefore elect to use the globular cluster luminos- 
ity functions of the two brightest ellipticals in Eridanus A, 
NGC 1400 and the central galaxy NGC 1407, to provide 
new estimates of their distances. In ^ and ^ we dis- 
cuss the observations and data reduction, while in §^ we 
present the GCLFs and distance determinations. In addi- 
tion, the inferred globular cluster radial profiles and scaled 
population sizes of these two galaxies are investigated in 
§^| and respectively, in search of any signature of the 
peculiar environment in which they are found. 

2. Observations 

The data for this study of the globular cluster systems 
(GCSs) of NGC 1400 and NGC 1407 were acquired as part 
of two separate observing runs. The first set of observa- 
tions was obtained on the nights of November 14 and 16, 
1993, with the Cerro Tololo Inter-American Observatory 
(CTIO) 4-m telescope. Use of the Tek 2048 x 2048 CCD 
camera at prime focus with an image scale of / .'46 per 
pixel yielded an image size of 15!7 x 15f 7. The readout 
noise was 6 electrons/pixel, and the gain was set to 3.2 
electrons per ADU. A total of three images, each one in- 
cluding both NGC 1400 and NGC 1407, were taken using 
the Washington Ti filter (Canterna 1976; Harris & Can- 
terna 1977). The seeing was mediocre at best, ranging 
from l'.'5 to 2'.'0 over the two nights. 

The second set of data was collected on January 2 and 

3, 1995, with the Canada-France-Hawaii 3.6-m telescope 
(CFHT). The Loral3 detector (2048 x 2048) was used at 
a nominal gain of 1.45 electrons/ ADU and with a readout 
noise of 9 electrons/pixel. This detector, mounted behind 
the re-imaging optics of the MOS instrument, yielded an 
image scale of / .'32 pix and a frame size of lOf 9 x 10!9. 
Six images, three each of NGC 1407 and NGC 1400, were 
obtained using the Kron-Cousins I-band filter. The seeing 
varied from l'.'l to l'/9 over the two nights. A log of the 
observations is presented in Table ^. In addition to the 
program fields, standard star fields, dome flatfields and 
bias frames were also obtained. 

3. Data Reduction 

Preprocessing of the raw frames included bias sub- 
traction and division by a mean dome flatfield expo- 
sure. Frames taken in comparable seeing were combined 
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Table 1 

Galaxies of Eridanus A with known velocities (Willmer et al. 1989). 



Galaxy 


R.A. 


Dec. 


V0 


B 




(1950) 


(1950) 


(km/s) 




NGC 1383 


3 35 23 


-18 30.1 


1948±19 


14.0 


NGC 1390 


3 35 37 


-19 10.3 


1215±33 


14.6 


NGC 1393 


3 36 23 


-18 35.4 


2185±26 


13.9 


NGC 1400 a 


3 37 16 


-18 51.0 


549±21 


12.1 


IC 343 


3 37 52 


-18 36.3 


1869±30 


14.5 


NGC 1407 b 


3 37 57 


-18 44.5 


1766±21 


11.1 


ESO 548-G48 


3 38 04 


-19 05.5 


1817±27 


14.4 


IC 346 


3 39 29 


-18 25.6 


1897±47 


14.1 


ESO 548-G79 


3 39 41 


-19 03.2 


2053±31 


14.4 


ESO 549-G02 


3 40 43 


-19 10.8 


1134±36 


15.0 



b The galaxies in this study 



Table 2 

Distance Measurements of NGC 1400 and NGC 1407 (Gould 1993). 



Method 


NGC 1400 


NGC 1407 


Reference 




Distance (Mpc) 


Distance (Mpc) 




SBF 


16.3 ± 1.0 


16.4 ± 1.0 


Tonry 1991 


D n - a 


(24 ± l)^ 1 


(19 ± lyr 1 


Faber et al. 1989 


X-ray/B 


(16±2)/i- 1 


(18 ± 1 


Donnelly et al. 1990 


(u - V)/V 


(40 ± 2) fr 1 


(30±2)/r 1 


Sandage & Visvanathan 1978; 



Visvanathan & Sandage 1977 
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Table 3 

Record of Observations. 



Frame 


Object(s) 


Date 


Start UT 


Exp. (sec) 


Filter 


Airmass 


Seeing (") 


TlA 


N1400/1407 


Nov. 14/93 


08:21 


600 


Ti 


1.534 


1.5 


Tib 


N1400/1407 


Nov. 16/93 


07:28 


900 


Ti 


1.294 


1.9 


Tic 


N1400/1407 


Nov. 16/93 


07:50 


900 


Ti 


1.394 


2.0 


Ia 


N1407 


Jan. 2/95 


08:48 


600 


/ 


1.447 


1.4 


Ib 


N1407 


Jan. 2/95 


09:01 


600 


/ 


1.289 


1.9 


Ic 


N1407 


Jan. 3/95 


08:01 


900 


/ 


1.345 


1.1 


Id 


N1400 


Jan. 2/95 


08:20 


600 


/ 


1.367 


1.5 


Ie 


N1400 


Jan. 2/95 


08:33 


600 


/ 


1.397 


1.5 


If 


N1400 


Jan. 3/95 


07:37 


900 


1 


1.290 


1.2 



to form composite images with improved signal-to-noise. 
The resultant composite images with substantially differ- 
ent seeing were reduced individually using the standard 
Sun/Unix implementation of IRAFf]. Figure |l| shows an 
example of a preprocessed frame in the T\ filter. 

The data reduction was largely automated with the aid 
of an IRAF script and the DAOPHOT subroutines (Stet- 
son 1987). A median filtering procedure was applied in 
order to create a model of the diffuse galaxy light which 
was then subtracted from the original image to better re- 
veal the underlying globular cluster system. The fitting of 
a model stellar profile (a "point spread function" - PSF) 
was carried out upon the filtered image to yield the un- 
calibrated photometry list. 

The fraction of objects successfully detected in each 
frame as a function of magnitude is known as the com- 
pleteness function (Harris 1990). Using the IRAF ADD- 
STAR task, a total of 10 3 artificial stars of known magni- 
tudes were added to each frame in 10 trials, thus increas- 
ing the population of objects in an individual frame by 
roughly 20% per trial. In order to compensate for the vari- 
ation in object detection with radial distance from the in- 
ner galaxy regions, separate completeness functions were 
determined for a series of annuli centred on the galaxy. 

The uncertainties associated with the computed com- 
pleteness function (/) are derived assuming a binomial 
distribution: 

2 „, /(I-/) m 
Of w — , (1) 

fladd 

where n a dd refers to the number of artificial stars added 




Fig. 1. — A preprocessed image of NGC 1400 (upper 
right) and NGC 1407 (lower left) in the Ti filter from 
CTIO. The size of the field is 15.'7 x 15.'7 with north at 
the bottom and east at the left of the frame. 



3 Image Reduction and Analysis Facility, distributed by the Na- 
tional Optical Astronomical Observatories, which is operated 
by AURA under contract with the NSF. 
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to a given bin (Bolte 1989). The detection limits (50% 
recovery levels) for the program frames are listed in Ta- 
ble |^. Using the known magnitudes of the artificial stars, 
we found that the systematic errors in the recovery process 
are negligible. An estimate of the internal uncertainty as 
derived from the scatter of the artificial stars was within 
0.01 <C a <. 0.05 for all frames over the magnitude range 
of interest. 

For galaxies as distant as NGC 1400 and NGC 1407, in- 
dividual globular clusters are indistinguishable from stars 
under typical seeing conditions. A galaxy's GCS reveals 
itself as a statistical overabundance of faint, star-like fea- 
tures concentrated around the galaxy. An effective means 
of identifying non-stellar objects employs the CLASSIFY 
routine (Harris 1990). This routine calculates radial mo- 
ments using weighted intensity sums around each object 
(see Harris et al. 1991 and references therein). We chose 
to use the effective radius parameter r_2 as a discrimina- 
tor: 



r-2 = 



g [j J -/(r 2 +0.5)] 

E 1 ! 



-1/2 



(2) 



where r is the radial distance of the j th pixel from the ob- 
ject centroid, and Ij is the intensity of the j th pixel after 
subtraction of the the local sky value. The summations 
were performed over those pixels within a maximum ra- 
dius of 10 pixels which had intensities greater than the 
3.5cr detection threshold. For objects of a given magni- 
tude, non-stellar objects such as faint background galax- 
ies tend to exhibit more extended wings and faint cores, 
thus deviating towards larger values of r_2. Upper limits 
were chosen such that roughly 95% of the artificial stars 
were correctly classified. In total, between 20% and 40% 
of the detected objects in each frame were identified as 
non-stellar and were removed from each photometry list. 

Photometric calibrations were performed in the stan- 
dard fashion. Fitting of transformation equations to the 
photometric magnitudes from short exposures of stan- 
dard star fields (Geisler 1996, PorteiQ) observed during 
the same nights as the program fields yielded: 

ii(ap) - Ti = 1.176 + 0.090 (C - Ti) + 0.048 X 



'(op) 



/ = 0.268 + 0.025 (V - /) + 0.114 X, 



where tu ap \ and i( ap ) are instrumental magnitudes, and Ti 
and I are the known magnitudes of the standards. Since 
only Ti and / frames were obtained for the program data, 
the colour indices C — T\ and V — I were set to 1.5 and 
1.0 respectively, typical values for globular clusters. The 
rms residuals for the transformations are ±0.05 in Ti and 
±0.04 in I. 



4. Globular Cluster Luminosity Functions 

4.1. Determining the GCLFs 

The observed globular cluster luminosity functions for 
NGC 1400 and NGC 1407 were obtained by binning the 
objects in the final photometry list into 0.4 mag intervals. 
We subtract the local sky LF (containing field objects and 
none of the GCS) from the observed LF (comprising both 
GCs and field objects) to reveal the GCLF for the target 
galaxy. No separate background frames were observed, 
hence the background density of objects was estimated 
from outer regions of the target galaxy frames where the 
GC density has dropped to zero (as determined from the 
radial profiles presented in the next section). 

The galaxy LF and that for the background region were 
divided by their respective completeness corrections in or- 
der to compensate for increasing detection incompleteness 
at fainter magnitudes. Subtraction of the corrected back- 
ground LF from the galaxy LF yields the globular cluster 
LF for each frame. The error bars on the luminosity dis- 
tributions reflect the uncertainty on the inferred number 
of objects per magnitude bin: 



n bs 



+ 



(1 - /)™obs 



Wadd/ 3 



(3) 



(Bolte 1989). Since the completeness (/) is the ratio of the 
observed number of objects (n bs) to the number inferred 
in a magnitude bin (n) , Eq. |l| allows us to rearrange Eq. ^ 
in terms of the known quantities for each bin to obtain: 



7l 1 + a/ 7 



(4) 



Once the final cluster luminosity distributions of each 
data frame were calculated, the observed GCLFs corre- 
sponding to the same galaxy and filter were averaged, 
weighted by the uncertainty on the number counts. Un- 
fortunately, for NGC 1400, the CTIO Ti frames in best 
seeing were not long exposures (see Table ||), while the 
deeper exposures were in poorer seeing than enjoyed at 
CFHT. For these reasons the NGC 1400 GCLF was not 
well delineated in Ti; we chose to work with the I ob- 
servations exclusively. The final background-subtracted, 
completeness-corrected and averaged GCLFs in T\ and I 
for NGC 1407 and in I for NGC 1400 are presented in 
Figures |, | and §. 

The globular cluster luminosity function is defined as 
the number of clusters per unit magnitude (m) and can 
be described to first order by a Gaussian function: 



(m) = Aexp 



-(m-m°) s 
2a 2 



(5) 



'Unpublished: IRAF PHOTCAL standards database. 



where A is a normalization factor, m° represents the 
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Table 4 

Completeness Limits. 



Frame 


50% Limit 


TlA 


7i = 22.6 


TlB/C 


T x = 22.4 


Ia 


I = 21.6 


Ic 


I = 22.2 


Id/e 


I = 22.4 


If 


I = 22.3 



1000 




1000 




Fig. 2.— The GCLF of NGC 1407 in Ti with Gaussian 
fits for a range of dispersions. 



Fig. 3.— The GCLF of NGC 1407 in 7 with Gaussian 
fits for a range of dispersions. 



turnover (peak) magnitude and a is the dispersion (Hanes 
1977). 

To fit the Gaussian function of Eq. ^, a nonlinear least- 
squares fitting procedure was applied. Since the data does 
not extend past the turnover magnitude, it was not possi- 
ble to leave all three parameters (A, m°, a) unconstrained 
during the fitting. Consequently, a series of values of 
1.0 < a < 1.5 were adopted and held fixed while the 
scale factor and turnover magnitude were permitted to 
vary. These values of a span the range typically found for 
GCLFs (Harris 1991). The best-fit results are shown in 
Tables ^, ^ and ^| Note that the uncertainties quoted here 
reflect the formal errors associated with the least-squares 
fit. 



4.2. Distance Determinations 

Comparing the observed brightness distributions for 
the globular cluster systems of NGC 1400 and NGC 1407 
with those known from local calibrators, we can calculate 
distances to these galaxies. With an adopted absolute 
turnover magnitude of = —7.27 ± 0.23 (a mean value 
for the 10 largest galaxies in Table 2 of Harris 1991) and 
assuming a = 1.2 for the GCLF shapes, we estimate dis- 
tances of 17.6 ±3.1 Mpc to NGC 1407 and 25. 4 ±7.0 Mpc 
to NGC 1400 (Table ||). Note that the turnover magni- 
tudes have been corrected for foreground absorption. For 
NGC 1407, we adopted galactic extinction corrections of 
Atj = 0.11 and Ai = 0.06, based on a value of A_g = 0.17 
(RC3). NGC 1400's value of A B = 0.14 yielded correc- 



6 



Table 5 

GCLF fit parameters for NGC 1407 in Ti. 



a 


A 


T° 


x 2 


1.0 


241 ± 15 


22.92 ±0.07 


5.2 


1.1 


264 ± 17 


23.13 ±0.07 


3.4 


1.2 


296 ± 22 


23.37 ±0.08 


2.8 


1.3 


336 ± 29 


23.64 ± 0.09 


3.0 


1.4 


389 ± 41 


23.93 ±0.11 


3.6 


1.5 


457 ± 59 


24.24 ±0.13 


4.5 



Table 6 

GCLF fit parameters for NGC 1407 in I. 



a 


A 


7° 


x 2 


1.0 


313 ± 33 


22.72 ±0.11 


6.3 


1.1 


358 ± 45 


22.97 ± 0.12 


4.9 


1.2 


418 ± 64 


23.25 ±0.15 


4.3 


1.3 


498 ± 92 


23.59 ±0.17 


4.3 


1.4 


615 ± 133 


23.89 ±0.20 


4.5 


1.5 


749 ± 194 


24.24 ± 0.23 


5.1 



Table 7 

GCLF fit parameters for NGC 1400 in I. 



a 


A 


1° 


x 2 


1.0 


111 ±25 


23.22 ±0.21 


15.4 


1.1 


130 ± 37 


23.50 ±0.25 


15.6 


1.2 


157 ±53 


23.80 ± 0.30 


15.9 


1.3 


194 ± 79 


24.13 ±0.35 


16.1 


1.4 


244 ± 118 


24.49 ± 0.41 


16.4 


1.5 


316 ± 178 


24.88 ± 0.47 


16.7 
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100 



-e- 



10 r 




19 



20 



21 



22 



23 



24 



Fig. 4.— The GCLF of NGC 1400 in I with Gaussian 
fits for a range of dispersions. 



tions of Atj = 0.09 and Ai = 0.06 for this galaxy. The 
distances we have calculated in Table [s] are inconsistent 
with the hypothesis that NGC 1400 is a foreground galaxy 
roughly 3 times closer than NGC 1407, as implied by its 
anomalously low recessional velocity. 

Further experimentation shows that the choice of dis- 
persion (a) has little bearing on the relative distances ob- 
tained for these galaxies. The distance ratios fall within 
the range d(iV1400)/d(7V1407) = 1.31 ±0.04 regardless of 
adopted a, suggesting that NGC 1400 is in the background 
with respect to NGC 1407. Adopting the worst-case sce- 
nario in which a = 1.5 for NGC 1407 and a = 1.0 for 
NGC 1400, we find that d(A1400)/d(iV1407) = 0.63, still 
placing NGC 1400 considerably further than the distance 
implied by its low recessional velocity. From this we may 
again conclude that, within the reasonable range of a, our 
observations of the GCS of NGC 1400 are not consistent 
with it being a factor of 3 closer than NGC 1407. 

If NGC 1400 were at a distance of 5.49/i." 1 Mpc, the 
turnover magnitude of the globular cluster luminosity 
function would be expected to appear at Tf 1 = 21.6 and 
1° = 21.0, assuming h = 0.8. With NGC 1400 GCS lu- 
minosity distribution data down to completeness limits of 
Ti = 22.6 and I = 22.4, our GCLFs should thus extend 
well beyond the turnover in such a case. Since the ob- 
served GCLFs for this galaxy do not show any evidence 
of having reached the turnover magnitude, this reaffirms 
our conclusion that NGC 1400 is not a foreground galaxy 
one third the distance of Eridanus A. 



5. Cluster Radial Profiles 

The surface density distribution of GCs plotted as a 
function of galactocentric radius reveals the projection of 
the spatial structure of the GCS. The annular bins used 
to construct the profiles were 75 pixels wide and concen- 
tric on the galaxy centroids. The number counts in each 
radial bin, corrected for completeness and divided by the 
total area of the annulus sampled, were taken to represent 
the surface density at the geometric mean annular radius 
(r — v / nnT'out). Subtracting the local background density 
(°bgd), found by adopting the mean value of the surface 
density beyond the radius at which the radial profile is 
no longer decaying, we obtain the distribution of the GC 
population alone. 

Non-linear least-squares fitting of the following rela- 
tions were performed on the radial profiles: 



logCTGC = a + b r 1/4 
log aac = c + d log r. 



(6a) 
(6b) 



The first model profile is a form of the empirical de Vau- 
couleurs law, and the second represents a scale-free power 
law (Harris 1986). The curves superimposed upon the ra- 
dial profiles shown in Figures ^ to ^ represent the best fits 
to the data with the parameters and formal errors listed 
in Table |^. In all of the radial profile plots, we can clearly 
identify the presence of a centrally concentrated globular 
cluster system. 

Surface intensity profiles of the galaxy halo light are 
also provided on the GC radial profile plots. These halo 
profiles were determined by fitting elliptical isophotes over 
the galaxy region, yielding the relative intensity of light as 
a function of radius. After subtracting the local sky level, 
the halo intensity was arbitrarily scaled to the amplitude 
of the GCS radial distribution for easy comparison. 

6. Specific Frequencies 

The specific frequency (Sn) of a galaxy is defined as the 
number of globular clusters per unit halo light, normalized 
to an absolute magnitude of My = 15: 



S N = N t x l °- 4 < M v+i5) 



(7) 



where the total population of GCs surrounding a galaxy is 
given by Nt (Harris & van den Bergh 1981). This param- 
eter conveniently reflects the size of the GC population, 
independent of galaxy luminosity. 

To calculate specific frequency we must first estimate 
the total number of clusters which make up the GCS of the 
galaxy. We use the radial density profiles from to count 
the total number of the clusters surrounding the galaxy 
down to the limiting magnitude (N b s ). We divide this 
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Table 8 

Galaxy distances assuming a — 1.2, corrected for absorption. 



Galaxy 






m 






V° 




(V°- 


-M°) 


d (Mpc) 


NGC 1407 


r o = 


= 23.26 


± 


0.08 


23.71 ± 


0.08 


30.98 


±0.31 


15.7 (til) 




7° = 


= 23.19 


± 


0.15 


24.19 ± 


0.15 


31.46 


±0.38 


19-5 (til) 


NGC 1400 


7° = 


= 23.75 


± 


0.30 


24.75 ± 


0.30 


32.02 


±0.53 


25.4 (tU) 



Table 9 
GC radial profile best-fit parameters. 







logaGC = 


a + br 1/4 


logCTGC = 


c ± d log r 


Galaxy 


Filter 


a 


b 


c 


d 


NGC 1407 


Ti 


4.70 ±0.56 


-3.01 ±0.47 


1.77 ±0.12 


-2.23 ±0.34 




I 


4.33 ±0.63 


-2.83 ±0.63 


1.46 ±0.11 


-1.74 ±0.38 


NGC 1400 


Ti 


4.62 ± 1.27 


-3.71 ± 1.22 


0.90 ±0.12 


-2.35 ±0.78 




I 


3.79 ±0.41 


-2.81 ± 0.46 


0.96 ±0.09 


-1.57 ±0.34 



result by the fraction of the GCLF which was observed 
to infer the total number of clusters in the system. To 
define the shapes of the GCLFs, we adopt a Gaussian form 
(Eq. |H|) assuming a — 1.2 and the best-fit parameters from 
Tables |, | andfg. 

Absolute V magnitudes are also required in order to 
determine Sn- According to RC3, NGC 1407 has a cor- 
rected total blue magnitude B T — 10.5 ± 0.2 and colour 
(B - V) T = 0.97 ± 0.01, yielding V T = 9.5 ± 0.2. For 
NGC 1400, the RC3 values are B T = 11.87 ± 0.13 and 
(B - V) T = 0.92 ± 0.01, giving Vi = 10.95 ± 0.13 for 
this galaxy. To translate V T into absolute magnitude M v 
we must adopt distances to the galaxies in question. For 
both galaxies, we assume a distance of 20.5 ± 1.2 Mpc 
based on the mean recessional velocity of the Eridanus A 
galaxies and Ho = 80 km s _1 Mpc -1 . The results of the 
specific frequency calculations appear in Table [lo| The 
specific frequencies calculated for NGC 1407 in the two fil- 
ters are in agreement within their respective uncertainties, 
and we therefore adopt a mean value of 5V = 4.0± 1.3 for 
NGC 1407 along with that of S N = 5.2±2.0 for NGC 1400. 

7. Discussion 



7.1. The Question of Distance 

The primary motivation of this study was to determine 
whether NGC 1400 is at a distance comparable to that of 
NGC 1407 and the rest of Eridanus A, or if it is instead 
a foreground galaxy as implied by its anomalously low 
recessional velocity. The globular cluster luminosity func- 
tions are not consistent with the foreground placement of 
NGC 1400, which indicates that this galaxy has a large 
peculiar velocity with respect to Eridanus A, the origin of 
which remains unknown. 

We may now add globular cluster systems to the evi- 
dence that NGC 1400 does indeed lie at roughly the dis- 
tance of Eridanus A, and should be counted as a member 
of the cluster. The distances we have calculated assuming 
a = 1.2 for the GCLF are slightly higher but still within 
two standard deviations from those derived by Tonry 
(1991) of 16.3 ± 1.0 Mpc (NGC 1400) and 16.4 ± 1.0 Mpc 
(NGC 1407), and are well within the range of potential 
distances previously quoted in the literature (see Table ^) . 

7.2. GCS Shapes and Specific Frequencies 

The spatial profiles of the globular cluster systems of 
NGC 1400 and NGC 1407 shown in §| seem rather unex- 
ceptional in shape. The radial distributions in Figures ^ 
and reveal that the projected density of NGC 1407's 
GCS falls off with the same slope as the halo light intensity 
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Table 10 
Specific Frequencies 



Galaxy Filter iV bs Nt Mv Sn 



NGC 1407 Ti 629 ± 76 2296 ± 279 -22.06 ± 0.33 3.4 ±1.1 
I 556 ±61 2985 ±327 -22.06 ± 0.33 4.5 ±1.4 
NGC 1400 I 106 ±32 922 ± 280 -20.61 ± 0.26 5.2 ± 2.0 



t i i i i r 



10 r 



1 r 




b = -3.01 
d= -2.23 
halo 



100 




(arcmin) 



r (arcmin) 



Fig. 5. — The background-subtracted radial density pro- 
file of the NGC 1407 GCS (in T x ) with de Vaucouleurs 
and scale-free power law fits provided (slopes given by b 
and d respectively). Also shown is the galaxy halo inten- 
sity profile. The background density level was found to 
be Obgd = 5.9 ± 0.3. 

for this central elliptical galaxy of Eridanus A. Figures 
and | indicate that the light profile for the NGC 1400 
halo may be slightly steeper than the cluster density dis- 
tribution, but this effect is not uncommon. It is somewhat 
unusual for GCS radial profiles to follow the halo light pro- 
file as closely as seen for NGC 1407; generally, the cluster 
surface density is more distended (Harris 1991). Accord- 
ing to Merritt (1983, 1984), the cluster population and 
halo structure of a system dominated by dark matter will 
remain largely unchanged (to first order) , not significantly 
altered by galaxy interaction processes. We are therefore 
likely seeing the cluster distributions as they were origi- 
nally formed, with GCSs which are roughly as centrally 



Fig. 6.— The NGC 1407 GCS radial density profile (in 
I), as in Figure ^. The background density was a hgd = 
9.8 ±0.6. 

concentrated as the halo stars. 

There exists an apparent relationship between the in- 
tegrated V magnitude of a galaxy and the shape of its 
GCS density profile (Harris 1986). The parameter gener- 
ally used to describe the shape of the radial distribution 
is the slope of the logarithmic profile. This parameter was 
given the label "d" in the power-law fitting of Section |B], 
but is commonly referred to as a. A plot of the slope 
parameter as a function of galaxy magnitude is shown in 
Figure ^| with data and references presented in Table [ll| 
A linear fit to the data yields the relation: 

a = -7.93 - (0.29 ± 0.03)My (8a) 

This is consistent with the relation as determined in Harris 
1993 (see correction appearing in Kaisler et al. 1996). 
The two galaxies in this study follow the overall 
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Table 11 

Correlation between radial density and galaxy luminosity. 



Galaxy 


Type 




a 


Source 


Milky Way 


Sbc 


-21.3 


-2.5 ±0.26 


H81, H76 


NGC 147 


dE 


-15.0 


-3.5 ±0.5 


H86, HR79 


NGC 185 


dE 


-15.2 


-3.5 ±0.5 


H86, HR79 


NGC 205 


dE 


-16.5 


-3.5 ±0.5 


H86, HR79 


NGC 224 


Sb 


-21.7 


-2.4 ±0.4 


H81, HR79 


NGC 524 


E/SO 


-22.1 


-1.7±0.12 


HH85 


NGC 1052 


E 


-20.8 


-2.26 ±0.27 


HH85 


NGC 1374 


E 


-19.8 


-1.8 ±0.3 


KKHRIQ96 


NGC 1379 


E 


-20.0 


-2.1 ±0.6 


HH86b, KKHRIQ96 


NGC 1387 


SO 


-20.2 


-2.2 ±0.3 


HH86b, KKHRIQ96 


NGC 1399 


E 


-21.2 


-1.54±0.15 


BHH91, KKHRIQ96 


NGC 1400 


E/S0 


-20.6 


-1.96 ±0.56 


This Paper 


NGC 1404 


E 


-20.8 


-2.2 ±0.5 


HH86b 


NGC 1407 


E 


-22.1 


-1.99 ±0.36 


This Paper 


NGC 1427 


E 


-20.0 


—2.0 ± 0.3 


KKHRIQ96 


NGC 3115 


SO 


-21.1 


-1.84 ±0.23 


HH86a 


NGC 3311 


cD 


-22.4 


-1.24 ±0.17 


H86 


NGC 3842 


E 


-23.2 


-1.2 ±0.4 


BH92 


NGC 4073 


cD 


-23.1 


-0.95 ±0.3 


BH94 


NGC 4278 


E 


-20.8 


-2.26 ±0.23 


Hv81 


NGC 4365 


E 


-21.6 


-1.15 ±0.25 


HAPv91 


NGC 4472 


E 


-22.9 


-0.89 ±0.13 


HAPv91 


NGC 4486 


E 


-22.7 


-1.55 ±0.08 


H86 


NGC 4565 


Sb 


-21.6 


-2.33 ±0.61 


FHPH95, H81, vH82 


NGC 4594 


Sa 


-22.1 


-2.0 ±0.2 


HHH84 


NGC 4649 


E 


-22.5 


-1.08 ±0.10 


HAPv91 


NGC 4944 


E 


-21.04 


-1.06 ±0.38 


FHPH95 


NGC 5018 


E 


-22.3 


-1.3 ±0.4 


HK96 


NGC 5128 


Ep 


-22.0 


-1.5 ±0.2 


HHHC84 


NGC 5170 


Sb 


-22.0 


-1.69 ±0.25 


FHHB90 


NGC 5481 


E 


-20.2 


-1.7 ±0.2 


MR95 


NGC 5813 


E 


-21.6 


-2.18 ±0.33 


Hv81 


NGC 7768 


cD 


-22.9 


-1.32 ±0.32 


HPM95 


UGC 9799 


cD 


-23.4 


-1.39 ±0.30 


HPM95 


UGC 9958 


cD 


-23.4 


-1.14 ±0.50 


HPM95 
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Fig. 7.— The NGC 1400 GCS radial density profile (in 
Ti), as in Figure ^. The background density was o"b g d = 
5.7 ±0.4. 



Fig. 8.— The NGC 1400 GCS radial density profile (in 
I), as in Figure ^. The background density was cr bg d = 
10.4 ±0.7. 



trend whereby higher-luminosity galaxies exhibit more ex- 
tended globular cluster systems. The measured slopes, 
a(NGC1407) = -1.99±0.36 and «(NGC1400) = -1.96± 
0.56, compare reasonably well with the predicted values 
of a = —1.53 and —1.95 respectively. If we do not in- 
clude the data points for the dwarf ellipticals or the spiral 
galaxies, a linear regression yields 

a = -7.55 - (0.27±0.05)Mv. (8b) 

The omission does not significantly alter the relation; 
both of the linear fits above are shown in Figure ^| 

The specific frequencies calculated for NGC 1407 (Sn = 
4.0 ±1.3) and NGC 1400 (S N = 5.2±2.0) are not unusual. 
Normal elliptical galaxies generally have Sn ~ 2 — 5, with 
variations from this range most likely attributable to dif- 
ferences in galaxy environment and GC formation mech- 
anisms. 

7.3. GC Colours 

Although the Ti and I filters do not provide a sufficient 
baseline to obtain cluster metallicities, we can estimate 
mean colours for the globulars surrounding NGC 1407 
and NGC 1400. Using the regions of overlap between 
the Ti and / images, we determine Ti — I colours for 
the GCs detected in both filters. The colour distributions 
are shown in Figure [To| the lack of a well-defined peak 
in the NGC 1400 GCS colour distribution is likely due 
to the poor quality of the Ti photometry for this galaxy, 
with its smaller cluster population. NGC 1407's cluster 



population is found to have a mean of Ti — I = 0.53 and 
a median value of 0.54, while for NGC 1400 we obtain 
mean and median colours of Ti — I = 0.54 and 0.55, re- 
spectively. These values are consistent with typical GC 
colours observed in other galaxies, a finding which con- 
firms the identification of the GCSs. 

8. Conclusions 

The principal results of this study of the globular clus- 
ter systems of NGC 1400 and NGC 1407 can be summed 
up in the following points: 

1. From the shapes of the globular cluster luminosity 
functions, we determine distances to these galaxies 
which place NGC 1400 at or beyond the distance 
of the Eridanus A group. This finding is in agree- 
ment with conclusions made using other methods 
previously cited in the literature. 

2. The shapes of the GCS radial density profiles and 
the specific frequencies of the two systems reveal 
no obvious abnormalities. This implies that, if Eri- 
danus A is as dominated by dark matter as its es- 
timated M/L value indicates, no anomalies are ev- 
ident from the GC spatial distributions and popu- 
lation sizes of its two largest galaxies. 

With a distance comparable to that of NGC 1407, 
NGC 1400 must have a high peculiar velocity in order 
to account for its exceptionally low radial velocity. Gould 
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Fig. 9. — The relationship between slope of the GCS 
radial distribution and total galaxy V magnitude. The 
solid line represents a linear fit to the data including all 
data points (Eq. [5a|), and the dashed line indicates the 
fit if one does not include dwarf or spiral galaxies in the 
sample (Eq. ^b|). The data points enclosed by large open 
squares highlight the values obtained for NGC 1400 and 
NGC 1407 in the present study. 

(1993) demonstrates that if the distance of NGC 1400 is 
consistent with that of Eridanus A, it must be bound to 
the sub-cluster due to a lack of other nearby mass concen- 
trations large enough to generate its high peculiar motion. 
The exact cause for the large peculiar velocity remains as 
yet unknown. It is possible that NGC 1400 has a large 
component of its velocity moving it towards the core of 
Eridanus A (ie: it has a large transverse velocity), reduc- 
ing its net radial velocity. Perhaps the velocity dispersion 
of the cluster has been severely underestimated since we 
only have velocity data for 10 of the 50 or so members, and 
that by some coincidence the other galaxies in the sample 
have significant velocity components perpendicular to our 
line of sight. If this is the case, NGC 1400's high peculiar 
velocity might not be particularly anomalous. 

The second enigma surrounding the Eridanus A sub- 
cluster is its abnormally high M/L ratio. It is possible 
that it is merely a dark cluster - many such clusters could 
exist which have so far avoided detection. More expansive 
surveys at higher limiting magnitudes in combination with 
reliable cluster-finding algorithms may reveal the presence 
of more dark clusters. The question remains: if indeed this 
cluster contains a great deal of dark matter, where did it 
come from, and why is NGC 1400 the only member (so 
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Fig. 10. — Colour distribution for the globular clusters 
surrounding NGC 1407 (solid line) and NGC 1400 (dashed 
line) . 

far) to show such a high peculiar velocity? Furthermore, 
why is there no evidence of strange effects on the radial 
distributions and population sizes of the member galaxy 
GCSs, for example, given this unusual environment? 

A more extensive analysis of the GCS of the two Eri- 
danus A galaxies NGC 1407 and NGC 1400 could be pro- 
vided by additional deeper multicolour photometry as well 
as spectroscopic observations. This data may contribute 
to a better understanding of the nature of the galaxies 
and their environment. There is an obvious lack of red- 
shift measurements for the majority of the Eridanus A 
galaxies (see Table 4 of Ferguson & Sandage (1990) for 
the complete membership list). This makes it difficult 
to determine an accurate M/L ratio for the cluster, as 
well as to derive any conclusions regarding the dynamical 
processes at work within Eridanus A. A more complete 
database of member galaxy velocities may shed some light 
on this dark cluster. 

The authors wish to extend their gratitude to the kind 
and helpful staff at CTIO and CFHT. This work was sup- 
ported by the Natural Sciences and Engineering Research 
Council of Canada in the form of operating grants (DAH 
and WEH) and a PGS A scholarship (KMP). 
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